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The study of fluctuations of particle multiplicities in relativistic heavy-ion reactions has drawn 
much attention in recent years, because they have been proposed as a probe for underlying dynamics 
and possible formation of quark-gluon plasma. Thus, it is of uttermost importance to describe the 
baseline of statistical fluctuations in the hadron gas phase in a correct way. We have performed 
a comprehensive study of multiplicity distributions in the full ideal hadron-resonance gas in differ- 
ent ensembles, namely grand-canonical, canonical and microcanonical, using two different methods: 
asymptotic expansions and full Monte Carlo simulations. The method based on asymptotic ex- 
pansion allows a quick numerical calculation of dispersions in the hadron gas with three conserved 
charges at primary hadron level, while the Monte-Carlo simulation is suitable to study the effect of 
. resonance decays. Even though mean multiplicities converge to the same values, major differences 

' in fluctuations for these ensembles persist in the thermodynamic limit, as pointed out in recent 

^ studies. We observe that this difference is ultimately related to the non-additivity of the variances 

^Sj . in the ensembles with exact conservation of extensive quantities. 

PACS numbers: 24.10.Pa, 24.60.Ky, 25.75.-q 



I. INTRODUCTION 

The statistical hadronization model has proven to be a very effective tool in describing average particle multiplicities 
^ID [ in high energy heavy ion reactions jlj] as well as in elementary particle reactions On the other hand, in this 
. model it would be possible to calculate also fluctuations of particle multiplicities once the status of the hadronizing 
' sources {clusters or fireballs) in terms of volume, mass, momentum and charges were known. Multiplicity and charge 
, fluctuations have been indeed proposed to be a good discriminating tool between quark-gluon plasma and hadron 
' gas 0, 0] provided that they survive the phase transition and the hadronic system freezes out in a non-equilibrium 
' situation. However, in order to properly assess the discriminating power of such observables, one should flrstly 
■ calculate fluctuations in a hadron gas by including all "trivial" effects, such as conservation laws, quantum statistics, 
1— H \ resonance decays, kinematical cuts etc. The effects of conservation laws on fluctuations in thermal ensembles have 
^ . been firstly addressed, in the perspective of heavy ion coUisions, in ref. 0- More recently, it has been pointed out 0,0 
■ ■ that in the canonical ensemble (CE) with exact conservation of charges, scaled second moment (scaled variance) of the 
multiplicity distribution of any particle does not converge to the corresponding GC value even in the thermodynamic 
limit, unlike the mean {3,12] ■ This was fairly understood among experts in statistical mechanics [Toj . but probably it 
has been shown explicitely for the canonical relativistic gas for the first time in refs. 0,0- Further deviations from the 
; ^ GC limit were found in the case of exact energy and energy-momentum conservation in the microcanonical ensemble 
d ^ (MGE) jr., J^j. Since in a heavy ion collision conservation of charges must be fulfilled, the difference between CE 
and GCE might have some impact on the estimated size of fluctuations in a statistical model. 

The calculations performed in these recent studies |H 001^3 were mainly concerned with simplified cases, such 
as pion gas and pion-nucleon gas. In the present work, we address the fiuctuations in the general multi-species hadron 
gas including all resonances up to 1.8 GeV mass and carrying three additive charges, that is baryon number _B, 
strangeness S and electric charge Q, in the CE and MCE. We give first a general formulation of the problem, and 
derive novel fluctuation formulae for the large volume limit in the CE and MCE. We also discuss the problem of the 
inequivalence between GCE and CE in the thermodynamic limit for scaled variance and we show that the ultimate 
reason thereof is the conceptual difference between additivity and extensivity |14|: while particle multiplicities are 
additive and extensive in both GCE and CE, variances are extensive (i.e. proportional to the volume) but they are 
non-additive in the CE, so that the scaled variance turns out to be a pseudo-intensive quantity (according to the 
definition proposed in ref. Q). 



X 



* Electronic address: becattini@fi.inf'n.it 
^Electronic address: antti. keranen@oul u.fi| 
tElectronic address: ferroni@fi.infn.it | 
§ Electronic address: ,gabbriellini@fi.infn.it| 



2 



The paper is organized as follows: in the next section we address the large volume limit of fluctuations and 
correlations in the canonical ensemble by the means of asymptotic expansions. Monte-Carlo calculations are described 
and compared to analytical ones in Sect. 3. In Sect. 4 we address the large volume limit of fluctuations in the 
microcanonical ensemble. In Sect. 5 we discuss the difference between thermodynamic limits of fluctuations in the 
different ensembles and relate them to fundamental properties of the scaled variance. In Sect. 6 we address the problem 
of measuring fluctuations of charged particles and calculate different proposed measures in the different ensembles. 
The results are summarized in Sect. 7. 



II. ASYMPTOTIC FLUCTUATIONS IN THE CANONICAL ENSEMBLE 

Following refs. [^Q) we describe fluctuations by means of the scaled variance of a multiplicity distribution: 



UJ ■ 



{N) 



(1) 



where N is meant to be the multiplicity of any hadron species, primary or final (i.e. after resonance decays) or the 
sum of an arbitrary number of hadron species (e.g. all negatively chraged). This is a finite quantity in the infinite 
volume limit because the difference between (A'^^) and (iV)^ depends linearly on the volume for large volumes. It 
is worth reminding that, if quantum statistics is neglected, the multiplicity distribution of any primary hadron is a 
Poisson, thus lo = 1. 

For the sake of simplicity, we will first keep our discussion at the level of the classical Maxwell-Boltzmann (MB) 
statistics. Indeed, none of our arguments is affected by this approximation, and at the end we will give the proper 
generalization to quantum statistics and discuss the corrections. In this framework, using the one-particle partition 
function: 
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and the fugacity Xj for each particle species j, the grand-canonical partition function can be written as: 



^Gc({A,})=nE 

and, consequently, multiplicities of different species are uncorrelated and Poissonianly distributed: 



(3) 
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Since the sum of random Poisson variables is still Poisson, this also holds for any given subset of particles, e.g. negative 
hadrons or baryons. 

In the canonical ensemble, the partition function does not factorize into one-species expressions because of the 
constraint of fixed charges. Let us consider a hadron gas with three abelian charges, i.e. baryon number B, strangeness 
S and electric charge Q. In the following, we will denote by Q = (Qi, Q2, Q3) = {B, S, Q) a vector with components 
these charges and by qj = (gij , (72,^, (73j) — {bj,Sj,qj) the vector of charges of the j'*'' hadron species I24|. The 
canonical partition function with charges Q can be written as: 
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where Wick-rotated fugacities Xj = exp[i^j Qi.j4'i\ ^-re introduced in the grand-canonical partition function Zqc- By 
setting Wi — exp[i(/)i], we may write Eq. (O as a triple integral over the unitary circle in the complex w plane: 
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The first and second moments of multiplicity distributions of a set h of hadron species can be calculated by inserting 
a suitable fictitious fugacity in the function Zqc, i.e. replacing Xj with XuXj in Eq. |SJ) if j S /i and taking the 
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derivatives with respect to A^, [if 
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Using these, the scaled variance can be written as the sum of a Poissonian term, i.e. 
term: 
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1, and a canonical correction 
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Therefore, in the canonical ensemble, the quantities appearing in the expressions of the moments of the multiplicity 
distributions are the canonical partition functions calculated for the difference between total charges and charges of 

The quantity within brackets in the equation above vanishes in the thermodynamic 

exp{—fl-Q/T) (see Eq. (fT^ below), jl being the vector of chemical potentials corresponding to 

the conserved charges Q. However, the factor {Nj)zi^(i) is proportional to V'^ and, if the difference between brackets 
has terms proportional to l/V, they could give a finite contribution to cj/j in the thermodynamic limit. In fact, we 
will show that if this limit is properly taken, this is generally the case. Let us first rewrite ■^q„^. by using Eq. @ as: 
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Pb = B/V, ps = S/V, Pq = Q/V being the baryon, strangeness and electric charge densities respectively. 

The thermodynamic limit of the moments will be calculated by using an asymptotic saddle-point expansion of the 
above integral in the large parameter V keeping pb,PStPq fixed. Note that the function / in Eq. H12|l will then be 
independent of the volume because so is according to Eq. On the other hand, g{w) does not depend on 

large parameters because hadron charges are limited to few units. Furthermore, with this choice, the function f{w) 
has a saddle point wq determined by the condition df{'w)/dwi — which is independent of the hadron charges qj in 
Eq. (|10() . It can be shown that the saddle point coincides with physical fugacities in the grand-canonical ensemble, 
i.e. Wq = (As, A5, Aq) T^. For instance: 
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and similarly for wb and ws- The last equation, along with the corresponding equations for B and S*, manifestly show 
that the solution wq must coincide with the fugacities in the GCE, as in this case Zk(i)W^g Wg" Wq = {Nk)GC- The 
leading factor in the saddle-point expansion of the canonical partition function Z^ will then be (see also Eq. l|Sn|) '): 
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i.e. proportional to the grand-canonical partition function itself. 

We will now proceed to calculate the large volume limits of the scaled variance lo in the pion gas and in the full 
hadron gas. The saddle-point asymptotic expansions of the canonical partition functions, for different set of charges, 
will not be derived in detail in the following. Their general expression is discussed and quoted in Appendix A. 



A. Illustrative example: classical pion gas 



The fluctuations in the canonical pion gas in Boltzmann statistics have been recently studied by Gorenstein et. 
al. Even though the pion gas (including 7r+ and 7r~, -k^ is irrelevant in this context) has only one conserved 

abelian charge, we apply the more general approach of saddle-point expansions, which is better suited in case of 
several conserved charges. This serves as a good introduction to the full hadron gas case. 

The partition function of the pion gas for a total charge Q can be obtained at the order 0{V~'^/'^) by using 
the general formula of the saddle-point expansion (|80|l quoted in Appendix A, with v — V , g{w) = l/w, /(w) — 
--pQ Inw + Zjr/V{w -f 1/w) and the saddle-point wq = Xq being the grand-canonical fugacity: 



^GC 
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where aiXq) and 7(Aq) are constants dependent on derivatives of the function / beyond second order in the saddle- 
point. Now /"(Aq) = \pQ + 2z^/V\q\ and defining C = Z^^/X'^"^'^ ^\/2ttV f'iXq) which is a constant factor 
independent of the pion charges, the previous equation turns into: 
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where we have used (tt )gc = ^tt/Xq, i.e. the expression of the grand-canonical mean multiplicity of negative pions, in 
the limit of Boltzmann statistics. Similarly, the asymptotic saddle-point expansion of the canonical partition function 
for a charge Q — Qi can be obtained by taking g{w) — Wq^ ^ in Eq. (|80|l : 
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Note that Eq. H16|) is, as expected, a special case of Eq. H17|l with qj = 0. From Eq. ([TJl and (|16|l . (|17|) we deduce that 
the mean multiplicity of pions converges to the grand-canonical value in the thermodynamic limit V ^ oo: 
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On the other hand, this is not the case for the scaled variance oj. In fact, the part in the parentheses in Eq. 
yields a term of order 0{V~^) stemming from the difference between the terms proportional to 1/{Q + 2(7r~)Gc) in 
Eqs. Hlt)|l . (11711 . On the other hand, the two terms proportional to a and 7(Aq), cancel out at the order 0{V^^) in 
the difference between Zq^q^^q^/Zq^q^ and Zq^q^/Zq. Specifically: 
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Inserting this and (|18|l into Eq. and writing Q — {tt^)gc — (tt )gc, we obtain the final expression of scaled 
variances for both net charge and total particle number distributions: 
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These nice and simple expressions reduce to the grand-canonical (Poissonian) limit 1 only in special cases, and never 
for all three quantities under the same conditions. From Eq. H2(J|I it can be seen that the Poisson value of 1 for the 
scaled variance is recovered in the thermodynamical limit for negatives if (7r+)GC S> (7r^)GC (and vice-versa), i.e. 
when the net charge of the system is largely positive and the oppositely charged particles become a small pollution; 
conversely, the scaled variance of the positive particles almost vanishes because (7r+)Gc — Qj which must not fluctuate. 



B. Full hadron gas 

We now turn to the general case of multi- hadron gas carrying total charges {Qi} — {B,S,Q). The asymptotic 
expansion of the canonical partition function involves diagonalization of the Hessian matrix of the function / defined 
in Eqs (p^l and (see Appendix A): 
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where qij is the I charge of the j hadron species. Because this is real and symmetric, the diagonalization is 
achieved through an orthogonal matrix A: 



H' = diag(/ii,/i2,/i3) = AHA^ 



(23) 



where hi are the Hessian eigenvalues. It should be emphasized that the hessian matrix H, as well as its eigenvalues and 
the diagonalizing matrix A, is not dependent on the volume, just like the function / in Eq. 1)12(1 . By using the general 
expression (|80|) quoted in Appendix A with / in Eq. (|12l) and g in Eq. it can be proved that the asymptotic 

expansions in the volume V of the partition function - reads: 
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SO that the usual grand-canonical expression is recovered for the average multiplicity in the large volume limit: 
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Furthermore, denoting: 



and using Eq. ((24(1 . it can be straightforwardly proved that: 
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It is not difficult to realize, from their definition, that the matrix M and the Cjk factors do not depend on the volume, 
i.e. they stay finite and non-vanishing in the thermodynamic limit. Therefore, we can recast Eq. ijSJ as: 



(28) 

oo, thus giving rise to 



where the second term on the right hand side of the above equation is finite in the limit V 
non trivial values of the scaled variance lo. 

In order to obtain a more intelligible expression of the extra-Poisson term in Eq. H28|l . let us consider the scaled 
variance ujj of a single hadron species in a simplified case, by neglecting baryon number and strangeness, i.e. consid- 
ering only electric charge. The hessian matrix is trivially diagonal and one is left with only one element in the matrix 



I.e. 
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Then, neglecting possible hadrons with two or more units of electric charge, we obtain the following expression for 
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where {h'^)Gc smd (/i~)gc ^I'e the mean multiplicities of positive and negative hadrons respectively. By plugging 
the last expression into Eq. (|28|l and taking into account that for a single hadron species j = k and = 1, we then 
obtain: 
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Likewise, it can be proved that, for the set of all negative, positive and charged hadrons, scaled variances read: 
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similarly to the classical pion gas case. The Eq. I|30() shows that the deviation from the Poisson statistics is proportional 
to the relative weight, in terms of multiplicity, of the species j with respect to all species carrying a non- vanishing 
value of the same charge. This nicely meets our physical intuition that in a system with many species, the distribution 
of single and rarely produced ones, is little affected by global conservation laws. Conversely, according to Eq. (|31|) . the 
scaled variance of the most inclusive sets are strongly affected by the conservation laws and can show large deviations 
from Poisson statistics. These general features hold in the more general case of three conserved charges, although 
the simple formulae (|3U|I and H31|l have to be modified. Nevertheless, they can be taken as "rules of thumb" to make 
estimates of scaled variances to a first approximation. 



C. Quantum statistics 



The generalization of previous expression to the case of quantum statistics is rather straightforward. The GC 
partition function reads: 
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where: 

with the upper sign for fermions and the lower for bosons. This yields the Hessian: 

1 
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The first and second moments become: 
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The middle term on the right hand side is the quantum correction to the Boltzmann statistics and is kept in GCE, 
whilst the rightmost term represents the CE correction to the GC expression. As far as this last term is concerned, 
it can be easily realized that the quantum hadron gas is equivalent to a classical hadron gas with an infinite number 
of particle species formed out of the actual ones by multiplying their mass and charges by a positive integer number. 
Thus, the same chain of arguments for the Boltzmann case can be followed in order to derive an asymptotic expression 
of uj: 



w/,. = 1 + 
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where Cjk has the same expression as in Eq. (I27|l . 



III. NUMERICAL RESULTS AND MONTE-CARLO CALCULATIONS IN THE CANONICAL 

ENSEMBLE 



In the previous section we have calculated the scaled variances of multiplicity distributions in the CE in the 
thermodynamic limit by using saddle point asymptotic expansion. We can now provide numerical values of these 
variances for several cases of interest for ultra-relativistic heavy ion collisions. In our calculation, all light-flavoured 
hadron species up to a mass ~ 1.8 GeV quoted in the 2002 issue of Particle Data Book 17] are included. The needed 
intensive input parameters for these calculations are the temperature T and the charge densities. The baryon density 
Pb is varied between and 0.3 fm^'^, while the strangeness density ps is set to zero and the electric charge density 
is set to pq = OApB, corresponding to the ratio Z/A of Pb-Pb and Au-Au collisions. The chemical potentials and 
fugacities are determined accordingly. 

The asymptotic scaled variances of the multiplicity distributions of primary baryons, strange and negative particles, 
as well as primary 7r~, p and K'^ for a completely neutral hadron gas are quoted in tabled The values in the CE in 
the Boltzmann approximation (B) are compared to the full quantum statistics ones in the CE (Q) and in the GCE 
(Q GC). It can be seen that lub, (^s Sind lo- are essentially 1/2 independently of the temperature, while cjtot — 1- 
This is not accidental, but a consequence of the vanishing of charge densities. To prove this, let us consider the scaled 
variance ojb+b of baryons and antibaryons. For this special case, it can be realized that the numerator on the right 
hand side of Eq. (|28|) . i.e. the extra-Poisson contribution to the scaled variance, vanishes. The reason is simply that 
the factors Cjk (see Eq. ()27f) ') enter in the sum with opposite values for a baryon-antibaryon pair with respect to a 
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TABLE I: Scaled variances for the full ideal hadron-resonance gas in the CE for T = 120, 160, 180 MeV and vanishing charge 
densities in the thermodynamic limit for different sets of primary hadrons (baryons, strange particles, negatives, all, pions, 
protons and kaons). In the statistics row, "B" stands for Boltzmann, "Q" for Quantum and "Q GC" for quantum grand- 
canonical. 



r [MeV] 


120 


160 


180 


statistics 


B 


Q 


Q GC 


B 


Q 


Q GC 


B 


Q 


Q GC 


U)B 


0.500 


0.500 


1.000 


0.500 


0.500 


1.000 


0.500 


0.500 


1.000 


ws 


0.504 


0.506 


1.004 


0.517 


0.520 


1.006 


0.524 


0.528 


1.007 


U>- 


0.502 


0.535 


1.066 


0.509 


0.536 


1.055 


0.512 


0.535 


1.045 




1 


1.061 


1.061 


1 


1.047 


1.047 


1 


1.038 


1.038 




0.603 


0.642 


1.088 


0.762 


0.824 


1.119 


0.831 


0.911 


1.131 


U)p 


0.892 


0.892 


1.000 


0.931 


0.931 


0.999 


0.942 


0.941 


0.998 




0.795 


0.804 


1.005 


0.865 


0.879 


1.013 


0.896 


0.913 


1.018 



baryon-baryon or antibaryon-antibaryon pair, while {Nj)Qc{Nk)GC has the same value for all above combinations if 
the gas is neutral. Since there is an equal number of baryon-antibaryon and baryon-baryon or antibaryon-antibaryon 
pairs, the sum in the numerator of Eq. (|28|l vanishes, thus lob+b = 1- In terms of the variance of the random variables 
Nb + Ng, we then have: 

{S{Nb + Nsf) = (Nb) + (Nb) = 2{Nb) (38) 

where we have used the neutrality condition. However, this condition holds for the random variables themselves and 
not just for their mean values, that is Nb = Nj^. Hence, the left hand side can be written as 4:{SNg), thence: 

{6{Nb+Nb?)=M6NI)=2{Nb) =^ ^B = J^ = l (39) 

i.e. the scaled variance of baryons is 1/2. The same argument could be carried over to ujs and uj± if there were 
only singly strange and charged hadrons. In fact, the slight difference from 1/2 for ujg seen in table ^ is due to 
multi-strange particles, while quantum statistics effects in primary pion production are the main responsible for the 
similar deviation of a;_ . 

For primary protons, pions and kaons, the scaled variance is not too far from its Poissonian limit; this is a con- 
sequence of the overall large number of species in the hadron gas, in accordance with the argument at the end of 
Subsect. lIIBI As pion and kaons and protons are the lightest charged, strange and baryon particles respectively, they 
are relatively more abundantly produced at low temperatures than at high temperatures and this is reflected in the 
growth of u! towards its Poisson limit going from T — 120 MeV to T = 180 MeV. 

TABLE IL Scaled variances for the full ideal hadron-resonance gas in the CE with quantum statistics for different sets of 
primary hadrons (baryons, antibaryons, strange, antistrange, positives, negatives, charged), in the thermodynamics limit, as a 
function of T and baryon density pB- Strangeness density ps ~ and electric charge density pq = 0.4p_B. 



r [MeV] 


120 


160 


180 


PB [fm-»] 


0.0 


0.15 


0.30 


0.0 


0.15 


0.30 


0.0 


0.15 


0.30 


LOB 


0.500 


0.005 


0.008 


0.500 


0.043 


0.015 


0.500 


0.195 


0.085 


Ug 


0.500 


1.000 


1.000 


0.500 


0.952 


0.980 


0.500 


0.799 


0.905 


LOS 


0.506 


0.474 


0.481 


0.520 


0.474 


0.465 


0.528 


0.507 


0.485 


LUg 


0.506 


0.493 


0.494 


0.520 


0.493 


0.486 


0.528 


0.514 


0.498 


UJ+ 


0.535 


0.255 


0.209 


0.536 


0.391 


0.316 


0.535 


0.464 


0.390 


LU^ 


0.535 


0.571 


0.566 


0.536 


0.555 


0.549 


0.535 


0.562 


0.558 




1.066 


0.649 


0.543 


1.055 


0.872 


0.749 


1.045 


0.974 


0.872 



In table ^ we quote the asymptotic scaled variances of particles and anti-particles for different temperatures and 
baryon densities, for full quantum statistics. It can be seen that lub and ujb are found to quickly converge to the limits 
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and 1 respectively with increasing ps, as predicted by the rule of thumb for the baryons analogous to Eq. I|31|l . 
On the other hand, the strange sector is dominated by the strangeness neutrality condition and is weakly dependent 
on baryon density. The fluctuations of positive hadrons are dominated by protons at high baryon fugacities. This 
is reflected in a monotonically decreasing value co^ with increasing ps- Conversely, on the basis of the rule of 
thumb H31() , one would expect an equally strong increase of uj- . Surprisingly, this is not the case because at high A b 
the A~ becomes the dominant negative hadron and that simple approximation fails because of the strong correlation 
between electric charge and baryon number. 



A. Monte-Carlo simulations 



The scaled variances determined by means of analytical calculations have been compared with those obtained 
through Monte-Carlo simulations. The basic idea of this method is to extract randomly ii'-uples {Nj} of multiplicities 
Nj for each hadron species j according to the multi-species multiplicity distribution of the canonical ensemble and 
averaging therafter. This method allows to determine numerically, with a finite statistical error, not only scaled 
variances but also higher order moments and, in general, to visualize the shape of the distributions. Furthermore, 
this method makes it possible to make calculations at final hadron level, taking into account resonance decays and 
thereby allowing a comparison of theoretical calculations with actual measurements. The multi-species multiplicity 
distribution in the canonical ensemble has been determined in the form of a cluster decomposition in ref. JJJ : 



pm}) = ^ 



HE n . 



(40) 



where {hn-"\ are partitions of the integers J^j in the multiplicity representations, i.e. such that Nj — '^nj=i 



Hj = Yln-=i hn^; and Zj(„^.) are given in Eq. 

In the limit of the Boltzmann statistics, the distribution H40|l reduces to a product of independent Poisson dis- 
tributions, one for each species, with the constraint of charges conservation S^sr w - • A direct sampling of the 

distribution H40(l is very difficult though. The most effective method is the importance sampling technique, in which 
each event (namely a if-uple {Nj}) is weighted by the ratio w of the true distribution P{{Nj}) and the actually 
sampled distribution R{{Nj}). The latter should be a distribution quickly and efficiently sampled and, moreover, as 
similar as possible to P{{Nj}) to minimize statistical errors. In our case, we have chosen R{{Nj}) as the product 
of unconstrained Poisson distributions, like in Eq. 10} . Their mean multiplicities are chosen to be those of the GCE, 
that is (Nj) = Zj(i)Aj, where the fugacities Xj are determined according to the saddle-point equations thereby, 
the mean values of the Poisson distributions in R{{Nj}) coincide with the actual CE average multiplicities in the 
thermodynamic limit. 

The average of any observable O after Ns random samples reads: 



Ns 
i=l 



(O) = witii = „ ; (41) 



^ Rm}^) 

1=1 



statistical errors can be calculated accordingly . 

Unfortunately, with this method, it is not possible to calculate observables straight in the thermodynamic limit 
because the simulation can be carried out only with a finite volume. Instead, one can study the variation of some 
observable of interest as a function of the volume, fixing total charges or charge densities, and estimate the thermo- 
dynamic limit by extrapolating. It should be pointed out that a too large volume cannot be used in order not to 
diminish too much the efficiency of the Monte-Carlo calculation. In fact, according to Eq. 141|) . only the if-uples with 
non-vanishing weight, i.e. fulfilling Q — Njqj actually contribute to the average. On the other hand, random 
samples are extracted from R{{Nj}), which is a product of Poisson distributions, and only a small fraction of them 
will meet the charge conservation constraint. This fraction is the efficiency of the sampling and reads: 

{Nj} 3 ^' 
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It can be easily proved, by developing the above expression, that the efficiency decreases as volume increases. In other 
words, the larger the volume, the more difficult is to catch a configuration which fulfills exactly the charge conservation 
constraint. Then, to keep the statistical error small, one should not increase the volume too much. Fortunately, the 
statistical error on most averages decreases as volume increases because of the increase in multiplicity of single events, 
so that a lower efficiency at larger volumes does not spoil the accuracy. 

We have run simulations of canonical systems with vanishing charges (neutral systems) and with charges equal to 
those of a pp collision, i.e. Q = 2,B = 2,S = (pp-like system) for T = 160 MeV and volumes V ranging from 5 to 
85 fm^ in steps of 10 fm^; for each point, 10^/y(fm'^) effective samples (i.e. fulfilling charge conservation) have been 
drawn from R{{Nj}) with efficiencies varying from 26% to 0.2% for neutral system and from 3% to 0.2% for pp-like 
system. In figs.^ 12 the obtained scaled variances are shown for negative and charged hadrons respectively, both at 
primary and final level. It can be seen that scaled variances slowly converge to their thermodynamic limits which are 
attained up to ~ 5% within V ^ 100 fm^ and that those of the finals are always larger than primaries'. The latter 
effect is mainly owing to neutral particles decaying into a pair of charged particles. 

We have also calculated the scaled variances in the canonical ensemble for larger systems and compared them with 
relevant analytical calculations. This is shown in fig. © where we have plotted u; for different sets of particles as a 
function of the baryon density ps, with pg — 0.4pB and ps = 0. Monte-Carlo results are shown as dots with errors 
bars and have been obtained by drawing 10^ effective samples with V = 200 fm'^ for baryon densities 0., 0.1, 0.2, 
0.3 fm^'^. Due to the large volume, which is needed to ensure the effective reaching of the thermodynamic limit, the 
efficiency of these runs is very low, in the range (1 — 6) • 10"''. It can be clearly seen that Monte-Carlo and analytical 
calculations are in excellent agreement. 



IV. ASYMPTOTIC FLUCTUATIONS IN THE MICROCANONICAL ENSEMBLE 



In principle, fluctuations in the microcanonical ensemble of the ideal hadron-resonance gas in the thermodynamic 
limit can be studied by applying the same saddle-point expansion described in Sect. 2 for the canonical partition 
function to the microcanonical partition function (MPF) , defined as '■ 



n = J2{hv\s*{p-Pop)s^^^jh^ 



(43) 



\hy) being localized multi-hadronic states within the volume V. However, for the sake of simplicity, we will confine 
ourselves to energy conservation only, i.e. momentum conservation will be disregarded. Furthermore, we will discuss 
only the case of Boltzmann statistics. Therefore, the operator (5'*(P — Pop) in Eq. (I43|l is henceforth understood to be 
replaced with S{E — i?op)- In this case, the partition function can be written as ,18j : 



ri/t — lim 



+ioo+e 

dwE 

-ioo-t-e 



dwf 



.0+ (27ri)4 

where the following convention is adopted: 



dws 4> dwQ ^ujg^ ^'^s^ ^^o^ '^exp'; 



(0) bj Sj qj 



in) 



(27r)3 



d^p e]w~'' 



, (44) 



(45) 



Sj being the energy of the j^^ particle. The integral expression of the microcanonical partition function 144|) looks very 
similar to the canonical partition function ©, with a further integration over the imaginary line we — ie. Insertion 
of the fictitious fugacities Xj allow to calculate the moments like in the CE in Eq. O. Denoting by: 



1 



nd. 



F 



(0) 



'Q (27ri)4 
it is quite straightforward to derive the moments: 



t;i(0) E-1 -C 
PjJ^E 



-^expjt/^P, 



(0) bi. Sz, Qk 



jeh Q 



j,k£h 



(46) 

(47) 
(48) 
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Consequently: 



-. = 1 + ^ '- (49) 



Q-qj 



jGh Q 

The whole procedure of saddle-point expansion for the canonical partition function can be carried over here for 
the microcanonical partition function. The saddle point equation in the variable we leads to the definition of the 
temperature T = 1//3 [3 

which is such that: 

E J ^ ' ^bj <Sj ^<^j _ 
y ^B^S - PE 

j 

T.i^.^^^'B^y^'Q-y i-W (50) 

The appearance of an exponential of a dimensional quantity in the equations will not cause any problem because will 
always have energies as exponents in the final expressions. The mean multiplicity (|47|) converges to the corresponding 
grand-canonical value in the thermodynamic limit (see Appendix B). Conversely, the term within brackets in the scaled 
variance expression H49|l is 0{V^^) yielding a non-trivial value in the thermodynamic limit. In fact, the following 
expression is obtained for lu for V ^ oo: 



E 



3 

e^ 



3 _ .... ._ --/3;^-i 3 

1 ^ j l,n=l 



[(li,kFj Ffc +qi,jFj Ffe J Moi + 2^ qijqm.kF^ A^^ A^^ A/,„ V . (51) 



where i^^"'' = F(")(e^) and, as usual, M;„ = J^^^q AuAm/hi, the index pertaining to the energy-related variable 
We- The details of the derivation of Eq. H51|) can be found in Appendix B. 



A. Numerical calculations 



The formula H51(l has been used to evaluate the scaled variance in the MCE in the thermodynamic limit. In 
table lllll we compare the scaled variance of several sets of hadrons with the corresponding value in the CE, for different 
temperatures. It is seen that that lu is further reduced owing to the additional constraint of energy conservation. This 
reduction is more considerable for larger sets of hadrons and attains its maximum for all primary hadrons {iVtot)- 

Also in the MCE, there are available Monte-Carlo algorithms to make numerical calculations of particle distributions 
both at primary and final levels. We have used an optimized importance sampling algorithm described in detail in 
ref. which is the state-of-art in this field. Unlike in the case of CE, Monte-Carlo simulations in the MCE of the 
full ideal hadron-resonance gas cannot be performed beyond values of ^ « 50 fm'^ and energies larger than « 20 GeV, 
with present computing power and numerical techniques . These values are not large enough to make an accurate 
numerical assessment of the thermodynamic limit. This can be seen from figs. 0)121 where the Monte-Carlo calculated 
scaled variances for increasing energies have been compared to their thermodynamic limit. Although at -E = 20 GeV 
central values sizeably differ from the thermodynamic limit, yet the overall trend is quite clear. It should be noted 
that, unlike in the asymptotic expansions, in the Monte-Carlo also momentum conservation was enforced, so that 
there might be a residual effect due to this difference in the calculation scheme. 

From a closer look at figs.01|Sl it can be realized that, unlike in the CE case, the scaled variance of negatives and 
charged is larger for primary than final hadrons especially at low energies and volumes. This might seem surprising, as 
our first intuition is that the decays of primary neutral particles into charged introduce a further source of fluctuation, 
thus leading to an increase of ui. This is indeed the case for the CE. However, in the MCE there is a further 
constraint, namely that the minimum number of primary particles is two. In a completely neutral MCE at very low 
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TABLE III: Scaled variances for the full ideal hadron-resonance gas in the MCE and CE for T = 120, 160, 180 MeV and vanishing 
charge densities in the thermodynamic limit for different sets of primary hadrons (baryons, strange particles, negatives, all, 
pions, protons and kaons). Momentum conservation and quantum statistics have been disregarded. 



T [Mo VI 
± livit; V J 


120 


160 


180 


dTi com nlo 


MCE 


CE 


MCE 


CE 


MCE 


CE 


LOB 


0.426 


0.500 


0.340 


0.500 


0.311 


0.500 


LOS 


0.357 


0.504 


0.349 


0.517 


0.346 


0.524 


LU^ 


0.307 


0.502 


0.310 


0.509 


0.304 


0.512 


f^ch 


0.610 


1 


0.603 


1 


0.584 


1 


Wtot 


0.257 


1 


0.207 


1 


0.172 


1 


LU^- 


0.523 


0.603 


0.736 


0.762 


0.815 


0.831 


UJ-p 


0.882 


0.892 


0.920 


0.931 


0.932 


0.942 




0.756 


0.795 


0.847 


0.865 


0.883 


0.896 



TABLE IV: Scaled variances for the full ideal hadron-resonance gas in the MCE for different sets of primary hadrons (baryons, 
antibaryons, strange, antistrange, positives, negatives, charged), in the thermodynamics limit, as a function of T and baryon 
density pB- Strangeness density ps = Q and electric charge density pq — OApB- Momentum conservation and quantum 
statistics have been disregarded. 



r [MeV] 


120 


160 


180 


PB [fm-^] 


0.0 


0.15 


0.30 


0.0 


0.15 


0.30 


0.0 


0.15 


0.30 


LOB 


0.426 


0.000 


0.000 


0.340 


0.036 


0.012 


0.311 


0.133 


0.064 


LOg 


0.426 


0.996 


0.998 


0.340 


0.808 


0.896 


0.311 


0.543 


0.686 


UJS 


0.357 


0.356 


0.362 


0.349 


0.341 


0.343 


0.346 


0.341 


0.340 


LUg 


0.357 


0.367 


0.369 


0.349 


0.360 


0.362 


0.346 


0.354 


0.358 


UJ+ 


0.307 


0.192 


0.172 


0.310 


0.246 


0.218 


0.304 


0.267 


0.242 


LU- 


0.307 


0.438 


0.464 


0.310 


0.374 


0.403 


0.304 


0.342 


0.369 


Wch 


0.610 


0.484 


0.440 


0.603 


0.562 


0.524 


0.584 


0.572 


0.551 


Wtot 


0.257 


0.141 


0.107 


0.207 


0.173 


0.148 


0.172 


0.163 


0.150 



energy particles are essentially produced only in pairs because the channels with 3, 4 etc. particles have neghgible 
probability |llj |. thus one can have zero (when two neutral particles are emitted) or one negative hadron (when 
one positive and one negative particle are emitted) at primary level. It is not difficult to realize that the relative 
probability of these two events is 1/2 if there is an equally large number of positive, neutral and negative species 
and that the corresponding scaled variance is 2/3, what seems indeed to be the limit of u- at small energy value in 
fig- El Conversely, at final level, if all neutral particles decayed into one positive and one negative particle, there would 
always be one negative particle and the scaled variance would vanish. The non-vanishing value is owing to neutral 
stable particles, neutral particle decaying into stable neutral particles (e.g. ir^ — > 77) and negative particles giving 
rise to more than one negative particle in their decays. Still, a strong reduction of scaled variance for final negatives 
is implied. A similar argument applies to charged particles. 



V. DISCUSSION 



We have found that the thermodynamic limit of the scaled variance is different in different ensembles. This effect 
has been understood for a long time in statistical mechanics. Though, from the previous derivations, the reader 
might have had the impression that this inequivalence between GCE, CE and MCE is a long-reaching consequence 
of a complicated analytical work not driven by a clear physical insight. Recently, it has been pointed out [T^ 
that variances are qualitatively different from particle multiplicities in that their proportionality to the volume is 
not "primordial" but arises from the difference of two quantities whose leading term is 0{V^) (see Eq. (Q). As 
a consequence, the behaviour of variance in the thermodynamic limit is determined by sub-leading terms in both 
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(A^^) and (A^)^ and different limits can be expected in different ensembles. We would like to point out here that the 
different behaviour in the thermodynamic limit can indeed be understood more simply and with more physical insight, 
by observing that, unlike particle multiplicities, variances are non-additive quantities in both the CE and MCE. 

The conceptual difference between extensivity and additivity has been recently discussed in ref. jl4| |. An additive 
quantity X is such that, if we split a general system into N subsystem: 

N 



hm y^X,^ Wxn Xce,mce{V) (55) 



On the other hand, an extensive quantity is such that the limit: 

lim — — X (53) 

N^oo N 

has a non-vanishing and finite value. If a quantity is additive, it is also extensive except for some exceptional case 
|l4j . Conversely, extensivity does not imply additivity. Quantities which are extensive but not additive are defined as 
pseudo-extensive and their corresponding limit x in Eq. (|53l) pseudo-intensive |l4j . It can be shown quite easily that 
additive quantities have the same thermodynamic limit in all ensembles. In fact, if we split a CE or MCE with a very 
large volume into a large number of N parts with volume V/N, each part is, by definition, a GCE with the rest of 
the system acting as a reservoir; this is just the way the GCE is introduced in most statistical mechanics textbooks. 
Consequently, any Xi, where i labels a subsystem in the CE or MCE, has the same value as in the GCE with volume 
V/N in the limit V, N oo. In other words, for any i: 

lim X,= lim ^^^^^^^ (54) 

V,N-*oo V.N^OD N 

If X is additive, then: 

N 

V,N^oo 

On the other hand, the left hand side of the previous expression also yields, according to 

lim J2 ^^^^ - lim XacEiV) (56) 

i=l 

The comparison between (|55f) and H5t)|) proves the equivalence between GCE and CE, MCE. Simplest examples of 
additive quantities are the energy and entropy for weakly interacting systems and particle mean multiplicities. 

The previous argument does not apply to the variance of particle multiplicity distribution. In fact, if we split a 
CE or MCE into N subsystems, the variance of any particle multiplicity distribution is not additive, as conservation 
laws involve non- vanishing correlations between different subsystems even for very large N. So: 

^VE^f (57) 

i=l 

The equal sign would apply only if the subsystems were completely independent of each other, which is the case only 
in the GCE. Thus, if variances are not additive, their GCE and CE thermodynamic limits may and, in general, will 
differ. Still, the variance itself is extensive because the limit H53|l yields a finite value in both CE and MCE, as the 
leading behaviour is 0{V). 

In conclusion, being non-additive, the variance is a pseudo-extensive quantity and the scaled variance is thus pseudo- 
intensive. For such quantities, the thermodynamic limit in the CE and MCE does not need to coincide with that in 
the GCE. 



VI. FLUCTUATIONS OF CHARGED PARTICLES 



Fluctuations of charged particle ratios on an event-by-event basis in heavy ion collisions have been suggested as 
probes of the prehadronic phase 0, Q and relevant measurements have been performed both at RHIC !l!^ and SPS 
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[20j . To start with, it is very important to stress that the mere variance of a ratio of extensive quantities is an ill-defined 
observable in statistical mechanics because it is not a (pseudo-) intensive quantity and vanishes in the thermodynamic 
limit simply because it is proportional to 1/V. Considering for instance the ratio Nj/Nk of particle multiplicities of 
two different species j and k we have: 




2cov(iV„7Vfc) (58) 



Since \cov{Nj, Nk)\ < y (SNj) (5N^) and (SN^) cx V, one is left with an expression which decreases at least pro- 
portionally to 1/V. Thus, in order to give a sensible fluctuation measure which does not vanish simply because the 
system gets larger, one should form some truly (pseudo-) intensive variable. 

Many such variables have been proposed to measure charge fluctuations in heavy ion collisions (for a review see ref. 
[21I), e.g. D [3, $Q m, and i/dyn Ej- Their definitions read: 

D -- 



t'dyn 




^Q = \lh^-V{Sq^) (59) 

/ 1 1 



where A^+, and A'ch is the number of positive, negative and charged particles respectively, Q ~ QjNj is the net 
charge and q is the charge of a produced particle. Note that in definition, in the first term the random variables 
are the numbers Nj of particles with a given charge, whereas the random variable in the second term is the charge 
of each particle itself. These variables are indeed related to each other and to the scaled variances of charged 
particles. It is not difficult to realize that the first two in Eq. H59|l are pseudo-intensive whilst the latter is not and 
should be multiplied by an extensive variable, e.g. (Nch) to make it such. Being pseudo-intensive, they have different 
thermodynamic limits in the GCE, CE and MCE. Therefore, much care is needed in comparing the measurements to 
the predictions of statistical mechanics because the effect of conservation laws is crucial in determining their values 
even for very large systems. 

However, this comparison is in general difficult because of additional source of fiuctuations which cannot be dis- 
regarded. Even if a statistical model framework would be essentially correct, there could be large fiuctuations of 
thermodynamic parameters (volume, temperature, baryon density etc.) from event to event, which, being super- 
imposed to the purely thermodynamic fluctuations, could swamp the thermodynamical fiuctuations. Furthermore, 
experimental measurements cover a limited kinematical window and this introduces a further complication. 

That said, it can be interesting to study the difference between the fluctuations of charged particles in the canonical 
and grand-canonical ensemble. In this respect, our Monte-Carlo method is especially suitable as it allows to predict the 
numerical values of the aforementioned variables taking into account all "trivial" effects including quantum statistics 
and resonance decays. Thus, we have calculated the values of these variables in the two ensembles at flnal hadron 
level for conditions relevant to heavy ion collisions at \/s^7v ~ GeV, i.e. T = 160 MeV, S = 0, Q/B = 0.4 and 
Pb = 0.2 fm^'^ by using Monte-Carlo simulations. In the CE, the volume chosen was 200 fm'^, which is large enough 
to ensure in practice the reaching of the thermodynamic limit. In order to show the effect of a limited kinematic 
acceptance, we have also implemented a toy dynamical model, giving each primary generated particle (according to 
thermal distributions) a random longitudinal boost in rapidity uniformly between —yh and yb- Though unrealistic, 
this model allows to understand the possible effect of measuring variables relevant to fluctuations over a finite rapidity 
window Ay. The calculated values of $q, D and (iVch)t'dyn are shown in figs. El [7||H1 respectively. It can be seen that 
all of them are strongly affected by the dynamical boost yb and the acceptance window Ay. Also, a relative strong 
difference is seen between CE and GCE. Yet, we note that, at least for <I>q and {Nch)i^dyn their CE value converge 
to the GC one for small rapidity acceptance. This is not a trivial feature because such a behaviour is expected if we 
select a subsystem in space and not in momentum space as we have actually done. In fact, this behaviour is not seen 
in D (see fig. 0). Altogether, we can conclude that the spread of these variables is considerable and making a fairly 
accurate estimate of the theoretical expectation for a hadron gas in chemical equilibrium requires at least taking into 
account exact charges (i.e. B, S, Q) conservation. 
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VII. SUMMARY AND CONCLUSIONS 



We have studied the thermodynamic limits of fluctuations of particle multiplicities in the full ideal hadron-resonance 
gas in the canonical and microcanonical ensemble taking into account the exact conservation of three charges (baryon 
number, strangeness, electric charge) and energy. The inequivalence between these ensembles and the grand-canonical 
ensemble in this respect is clearly seen in our calculations, which have been carried out by means of two independent 
methods: asymptotic expansions of the partition functions, giving rise to analytical formulae for scaled variances in 
the thermodynamic limit, and a full Monte-Carlo simulation. Excellent numerical agreement has been found between 
the two. 

The difference between the thermodynamic limits of the scaled variance in the MCE, CE and GCE ensembles has 
its roots in the non-additivity of the variance in the former two. In fact, we have shown that all additive quantities, 
like particle multiplicities, must have the same thermodynamic limit. 

We have finally discussed and calculated some variables used to measure charged particle fluctuations in the canon- 
ical and grand-canonical ensemble taking into account conservation laws and the decay of all resonances. It has been 
shown that their values is, in general, strongly affected by acceptance window and superimposed dynamical effects. 



APPENDIX A - SADDLE-POINT EXPANSIONS 



Generally, we want to find an asymptotic expansion of the complex c?-dimensional integral: 

d 



w- 

J-i 27ri 



.fc=i 



g{w) e 



i//(tu) 



(60) 



for large values of the real parameter v; w is a d-dimensional vector and (?(w) a smoothly varying function along the 
paths of integration Ffe. As it is well known, with large values of v, the dominant contribution to the integral comes 
from the small part of the path in the neighbourhood of the saddle point wq defined by 



Wo : 



df{w) 



dwk 







fc = 1, 



In this neighbourhood, one may expand /(w) in Taylor series up to second order as: 



1 9^/ 
f{w) ~ f{wo) + ^^{wi- woi) [wk - wok) 



i.k 



dwidwk 



(61) 



(62) 



According to the saddle-point or steepest-descent method [23, we choose a real integration variable tk for the path 
which is related to the original w through: 



Wk - wok = e"'"'tk 



(63) 



(j>k being a phase. If WQk lies on the real axis and the eigenvalues of the Hessian matrix are positive, as it is the case 
for canonical partition function, the steepest descent method prescribes that this phase ought to be 7r/2 if the so: 



Wk - wok = itk 



(64) 



This transformation corresponds to deform the original path into the line Wk — wok in the complex plane. Since 
only a small segment around the saddle point contributes to the integral value, we can extend the integration in tk 
from —00 to -l-oo leaving it essentially unchanged. By plugging the transformation 164|l into Eq. (|60|l and taking the 
expansion (|62|l for /, we can write: 



n 

U-1 



+00 



d<A 



9 



\ut^ Ht 



(65) 



where H is its Hessian matrix of the function /(w) for w — wo- The smooth function g{'w) is expanded into a Taylor 
series around the saddle-point: 



^ ^ yia a rik 
n=0 n llfc=l""'fe 



n 



Wn k — 1 



rikl 



(66) 
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where n = (ni, . . . , Ud) labels the sets of d integers such that n = J^k Thus Eq. (|65|l can be rewritten as: 



liu) ^ e 



lyf {wo] 



1 •„ d"9iw) 



d 



.fe=l 



+ 00 



dtk tl" 



e 2 



vt^ Ht 



So far, the problem is reduced to solving the integrals: 



n- 



.fe=l 



fife! J- 



e 2 



i^t^ Ht 



(67) 



(68) 



which is more easily achieved by diagonalizing the symmetric Hessian matrix in wq, where it is real. Thus, there must 
be an orthogonal matrix A such that: 

H' = diag(/ii,...,M =AHA^, (69) 
where hk are the Hessian's eigenvalues. By defining r = A^, the integral (|68l) can be rewritten as: 



, "fe. J -00 \ 1 



^UT H r 



.fc=l 

Now we can use the well known expansion: 

/ d Y'" 

( ^ ^ AmkTm ] 
\m=l / 



n / dr J ^ A™,r, 
.fe=i ^k- J-00 \^^^-^ 



(70) 



p'' p'' 



p'=(nfc) 



where p*''('^fe) = (Pi: • ■ ■ iPd) ^ '^^t integers such that X]-m=iPm = "fc- Therefore, the integral lf7n|l turns to 



(71) 



- n E 



ATk e-'^'^^'-^i n 



<L_ ( _A_ A 'Pm \ r + OO 

^mk 



E n n.fe, 



m 
m— 1 



(72) 



and, by rescaling the variables Tk^/hkV —>■ Tk'. 



yd/2 



E n 



p^(ni),...,p^(na}k=l \m=l ™ / V 



+00 



drtr: 



(73) 



Taking into account that X^mPm ~ and X]fc "-fc ~ ^'^d using known formulae for the rightmost gaussian integrals 
above, we are led to the following final expression of the integral (|67|) : 



liiy) e' 



(2,,).detH^^ 4rnti5<'= 



n E 



(pi + ...+p^-l)!! / A 



wo k=lp''(nk) ^k 



(Pfe + -+P^)/2 



n 



V,m=l ^™ 



(74) 



where, by convention, (—1)!! = 1 and iV!! = if TV is even. The Eq. H74|l is admittedly rather cumbersome, but it 
makes it clear that high order terms depend on increasing powers of l/i^, which are negligible for large values of v. 

Yet, the expansion H74|) is not a full asymptotic series for because the function / was expanded only at the 
second order. In order to estimate the corrections due to higher order terms in / expansion, one defines a function i] 
such that: 



i,k 



dw^dwk 



Obviously: 



f]{w) ^ O {{w - wof) 



(75) 



(76) 
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The series (|74() would turn into an exact asymptotic series provided that g is replaced with G(w) = g{w) ex.\i[vrj[wy\. 
However, this replacement affects the v dependence of the various terms in the series, because G{'w) depends now on 
V. This is most easily seen by expanding exp[vrj{'w)], i.e.: 

G{w) = g{w) ex\)[urj{w)\ ~ g{w) + i'g{w)'q[w) + —iPg[w)7f'{w) + -^v^ g{w)Tf {w) + . . . 

Because of (|76f) . the first two terms in (|74f) . with n = and n = 1, involving derivatives of G of order and 2, are 
unchanged, while the term n = 2, involving derivatives of G of order 4, is not: 



G{wq) ^ g{wo) 



d^G{w) 



d'^G{w) 



d^g{w) 



d^g{w) 



+ vgiwa) — 



dg{w) 



. dw, 



Ui-l 



(77) 



where we have used Eq. I|75(l . It can be seen that the rightmost two terms in the last equation are proportional to v, 
hence they introduce in the term with n = 2 of the series in Eq. I|74|) quantities of the order 0{v^^). Furthermore, 
even in the term n = 3 of the series a term in the sixth derivatives of G appears which turns out to be of the order 
0{v^^), that is: 



1 2 i^ .d^m) 

-v glwo) — ^ 

We can finally write the full expansion of /(j/) at order 0(1/"^) as: 

d 



(78) 



d^giw) 



„,i„i=4 rife^i 



n E 

wo k=l p''(nfc) 



2 dwkdw„ 

k,m—\ 



n 



k I 



E E^ 

n,|n|=4 i=l 



dg{w) 



dw. 



Wo 



1 d^fjw) 

2 r..M=6Uk=l9Wk 



d 



n E 

Wo ^"^1 p^irik) 



m—1 



TT ^mk 



vm=l ^™ 



n i: 

too '^=1 p''(nk) 



Pk 



i-m 



■J~j' ^mk 



ipl+...+pt)/2 \ 11 pk\ 

fl^ \m—l '"' 



(79) 



where, by convention |n| = n^. The above expression can be shortly summarized as: 

/ d , , \ d 



^yf{wo) 



1 



(27ri/)'*detH 



9{wq) + 



1 



IE 



d^g{w) 



2 /-^ dwkdwr, 

k.m—l 



E 



E"^ 



dg{w) 



dwi 



with 7 and a being constants dependent only on the function / and its derivatives. 



(80) 



APPENDIX B - SADDLE-POINT EXPANSION OF THE MICROCANONICAL PARTITION FUNCTION 

To work out the asymptotic expansion of the MPF in Eq. H44|l we define our / and g functions like in Appendix A 



as: 



9\w) 



We Wb wg w^ 



f{w) = pE In WE - Pb In wb ~ ps In ws - pq In wg + Pf'' {we)w^ w^ w^j. 



(81) 
(82) 
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Thus: 



lim 



1 



0+ (27ri)^ 



+ioo+€ 



dwE ^ dwB ^ dws ^ dwQ g{w) exp[V f{w)] 



(83) 



By using the expansion H8Q(I at the leading order for the functions 57^ and ^^i' defined in Eq. H46|l . plugging in the 

saddle point equations ()50|l defining temperature and fugacities, and proceeding like in Subsect. Ill Bl it is rather 
straightforward to obtain, in the limit of Boltzmann statistics: 



hm (Nh) = lim T/ V ■ ^ "^^ 



jeh Q 



(84) 



with Xj = A^^ Xq , that is the thermodynamic limit of mean multiplicities in the MCE is the corresponding value 
in the GCE. 



Let us now turn to scaled variances. Denoting fJ"^ = F)"^(e'3), with F)"\we) given by Eq. H45|l 



~ (27r)3 

and wq = (e^, As, As, Aq), the Hessian of / reads: 



Hoaiwo) 
Hoi{wo) 

Hiniwo) 



VX. 
1 



— e "l^qijFj Xj 



VXq,Xq„ 



QlSln + ^ qtj {qn.j - Sin) (^j)gC 



If we now take the general expression of the asymptotic expansion in 1)80(1 , it can be seen that the function g enters 
with both first and second order derivatives. However, since we are interested in calculating scaled variances, only the 
term involving second derivatives yields a non-vanishing contribution, just like in the canonical case. We will then 
confine ourselves to calculate the latter, taking into account that g varies according to which fi-related quantity we 
are interested in. In fact: 



Q-qj 



g{w) = Wj^^Wg^Wg^WQ^ 

^\ ]ri(0) -1 fcj-1 Sj-1 

g^-"(w) — I' J 

p(0) p(0) -1 6j+6fc-l Sj+Sfc-1 flj+gt-1 



-qj-ik 



gW^w) 



w, 



Q 



Nowe we can calculate the second derivatives of the functions g. Taking into account that: 



dwi 



iw) 



Wl 



it can be easily shown that 



(85) 
(86) 
(87) 

(88) 



for Eq. (|55|l. 



d^g{w) 



dwidwn 



{1 + Sir,){wQ)YH^oV lli^oV 



m=0 



/3 x"^^^ ( 2Ff + + Pf 



e ''A^ A^' A 



-2/3 



W( 



w 



dwEdwQi 
d^g^i\w) 



'0 



iqu~l)(Fr+Ff)e-% 



At/ Ag ^ (q,_^. _ 1) (g„ 



Sin) Ff^ Xq] XqI 



(89) 



(90) 
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for Eq. (|HSl, and 



dwEdwQ, 



dwQ^dwQ^ 



3 ^ k 



j k 



(91) 



ybj+bk-l -^Sj+Sk-l ^qj+qh-i 



+ - 1) (^r^r +Ff^Fn e-^A- 



for Eq. (|g7t. 

We can now use the general expression (|8()|l obtained in Appendix A to express the factor in parentheses in Eq. H49|l . 
Just Uke in the canonical ensemble case, the terms proportional to g and to the first derivatives of g within brackets 
in Eq. (|80|l with coefficients a and 7, do not contribute to the scaled variance in the limit V —> 00. Indeed, defining: 



4=0 

we obtain, at the order 0{V~^): 

Q-qk-qj Q-qk ^ ±_ ^ 



and 



d^g 



dwidwm. 



91: 



(92) 



Q-qj 



Q Ln=l 



g{wo) g('=)(wo) 



(93) 



Using Eqs. (|89|l . (|90() and (|91|l to work out the expression within brackets in Eq. (|93|l we finally obtain the scaled 
variance in the thermodynamic limit of MCE quoted in Eq. H51|l : 



lim LOh — I — ^ TTT-, 



k£h 



F^'^Ti'^ 

^(5^e Moo 



Y: [<l,kFf 'F':'+q,,Ffn'] -^M,n + E %3lm^kF':'^Q^^QlM,. 

1=1 Fa l,n=l 



(94) 



ACKNOWLEDGMENTS 



We gratefully acknowledge useful discussions with M. Gazdzicki, M. Gorenstein, K. Rajagopal, M. Stephanov, H. 
Touchette. 



[1] J. Cleymans and H. Satz, Z. Phys. C 57, 135 (1993); G. D. Yen, M. I. Gorenstein, W. Greiner and S. N. Yang, Phys. 
Rev. C 56, 2210 (1997); F. Becattini, M. Gazdzicki and J. Sollfrank, Eur. Phys. J. C 5, 143 (1998); P. Braun-Munzinger, 
D. Magestro, K. Redlich and J. Stachel, Phys. Lett. B 518, 41 (2001); F. Becattini, M. Gazdzicki, A. Keranen, J. Manninen 
and R. Stock, Phys. Rev. C 69, 024905 (2004); P. Braun-Munzinger, K. Redlich and J. Stachel, nucl-th/0304013 

[2] F. Becattini, Z. Phys. C 69, 485 (1996); F. Becattini and U. Heinz, Z. Phys. 76, 269 (1997); F. Becattini and G. Passaleva, 
Eur. Phys. J. C 23, 551 (2002) 

[3] S. Jeon and V. Koch, Phys. Rev. Lett. 85, 2076 (2000) 

[4] M. Asakawa, U. W. Heinz and B. MuUer, Phys. Rev. Lett. 85, 2072 (2000) 

[5] M. Stephanov, K. Rajagopal and E. V. Shuryak, Phys. Rev. D 60, 114028 (1999) 

[6] V. V. Begun, M. Gazdzicki, M. I. Gorenstein and O. S. Zozulya, Phys. Rev. C 70, 034901 (2004) 

[7] V. V. Begun, M. I. Gorenstein and O. S. Zozulya, nucl-th/0411003 

[8] J. Cleymans, M. Marais and E. Suhonen, Phys. Rev. C 56, 2747 (1997) 

[9] A. Keranen and F. Becattini, Phys. Rev. C 65, 044901 (2002) 



20 



H. Touchette, private communication. 

F. Becattini and L. Ferroni, Eur. Phys. J. C 38, 225 (2004) 

V. V. Begun, M. I. Gorenstein, A. P. Kostyuk and O. S. Zozulya, |nucl-th/0410044[ A. Keranen, F. Becattini, V. V. Begun, 

M. I. Gorenstein and O. S. Zozulya, J. Phys. G 31, S1095 (2005) 

J. Cleymans, K. Redlich and L. Turko, hep-th/0503133 

H. Touchette, Physica A 305, 84 (2002). 

F. Becattini and U. Heinz, Z. Phys. C 76, 269 (1997) 

F. Becattini, M. Gazdzicki, A. Keranen, J. Manninen and R. Stock, Phys. Rev. C 69, 024905 (2004) 
The Review of Particle Physics Ed. 2002, K. Hagiwara et al, Phys. Rev. D 66, 010001 (2002) 
F. Becattini and L. Ferroni, Eur. Phys. J. C 35, 243 (2004) 
J. Adams et ai, STAR CoU., Phys. Rev. C 68 044905 (2003) 

M. Rybczynski et al, NA49 Coll., J. Phys. Conf. Ser. 5, 74 (2005) and references therein; C. Ah et al, NA49 CoU., Phys. 
Rev. C 70, 064903 (2004) 

C. Pruneau, S. Gavin and S. Voloshin, Phys. Rev. C 66, 044904 (2002) 
S. Mrowczynski, Phys. Rev. C 66, 024904 (2002) and references th erein 
http:/ /www. math. ohio-state.edu~gerlach/math/BVtypset/nodell9. html" 

For the sake of clarity, it is worth stressing the difference between qtj, which is the ith charge of the hadron species j and 
Qj, which stands for its electric charge. Likewise, whilst Qi stands for the net ith charge of the system, Q is its net electric 
charge throughout the paper 



21 



1.3 
1.2 
1.1 
1 

0.9 
0.8 
0.7 ■ 
0.6 ■ 



0.5 F 

0.4 



Final ■ — b- 
Primary ■ — 
^ Asymptotic Primary 



m m 



10 20 30 40 50 60 70 80 90 
V [W] 



1.2 



0.8 



0.6 



0.4 



Final 
Primary 
Asymptotic Primary 



BBS 



m m m 



10 20 30 40 50 60 70 80 90 
V [fm^] 



FIG. 1: Scaled variances in the canonical ensemble of the full ideal hadron-resonance gas for negative hadrons at T = 160 MeV 
as functions of volume calculated with Monte-Carlo siunilations. Closed squares are for primary particles and open squares 
are for final particles after strong and electromagnetic decays. The horizontal solid lines indicate the thermodynamic limit 
y — > c» calculated with asymptotic expansions. Left panel: e"'"e~ charge configuration (Q=S=B=0). Right panel: pp charge 
configuration (Q=B=2, S=0). 
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FIG. 2: Scaled variances in the canonical ensemble of the full ideal hadron-resonance gas for charged hadrons at T = 160 MeV 
as functions of volume calculated with Monte-Carlo simulations. Closed squares are for primary particles and open squares 
are for final particles after strong and electromagnetic decays. The horizontal solid lines indicate the thermodynamic limit 
V —>■ oo calculated with asymptotic expansions. Left panel: e"'"e~ charge configuration (Q=S=B=0). Right panel: pp charge 
configuration (Q=B=2, S=0). 
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FIG. 3: Scaled variances in the canonical ensemble of the full ideal hadron-resonance gas for different sets of hadrons at T = 160 
MeV, S — and Q/B — 0.4 as functions of baryon density pB- Closed dots indicate the calculated values with Monte-Carlo 
simulations at primary level, open dots at final level with V — 200 fm'^. The lines depict the thermodynamic limits indendently 
calculated with the asymptotic expansion formulae. Left panel: charged, positive and negative hadrons; the arrows show the 
change in W- and u>ch from primary to final level. Right panel: baryons, antibaryons, strange and antistrange hadrons. 





FIG. 4: Scaled variances in the microcanonical ensemble of the full ideal hadron-resonance gas for negative hadrons at an 
energy density of 0.4 GeV/fm^ as functions of energy calculated with Monte-Carlo simulations. Closed squares are for primary 
particles and open squares are for final particles after strong and electromagnetic decays. The horizontal solid lines indicate 
the corresponding thermodynamic limit V ~> oo calculated with asymptotic expansions, the resulting temperature being ~ 160 
MeV. Left panel: e^e^ charge configuration (Q=S=B=0). Right panel: pp charge configuration (Q=B=2, S=0). 
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FIG. 5: Scaled variances in the microcanonical ensemble of the full ideal hadron-resonance gas for negative hadrons at an 

energy density of 0.4 GcV/fm"* as functions of energy calculated with Monte-Carlo simulations. Closed squares arc for primary 
particles and open squares are for final particles after strong and electromagnetic decays. The horizontal solid lines indicate 
the corresponding thermodynamic limit V ^ oo calculated with asymptotic expansions, the resulting temperature being ~ 160 
MeV. Left panel: e~^e~ charge configuration (Q=S=B=0). Right panel: pp charge configuration (Q=B=2, S=0). 




FIG. 6: Calculated $q (see text for definition) in the canonical and grand-canonical ensembles of the full ideal hadron-resonance 

gas at T = 160 McV and pB = 0.2 frri^"^, in the thermodynamic limit, at final hadron level for different random boosts of 
primary particles yt as a function of the acceptance rapidity window Aj/. 
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FIG. 7: Calculated D (see text for definition) in the canonical and grand-canonical ensembles of the full ideal hadron-resonance 

gas at T = 160 McV and pB = 0.2 fm"'', in the thermodynamic limit, at final hadron level for different random boosts of 
primary particles yb as a function of the acceptance rapidity window Aj/. 




FIG. 8: Calculated (A'ch)j'dyn (see text for definition) in the canonical and grand-canonical ensembles of the full ideal hadron- 
rcsonanco gas at T = 160 McV and pB = 0.2 fm^^, in the thermodynamic limit, at final hadron level for different random 
boosts of primary particles yt as a function of the acceptance rapidity window Aj/. 



